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Abstract 

A parameter optimization framework is presented to solve the problem of partitioning a central- 
ized controller into a decentralized hierarchical structure suitable for integrated flight/propulsion 
control implementation. The controller partitioning problem is briefly discussed and a cost function 
to be minimized is formulated, such that the resulting “optimal ,r partitioned subsystem controllers 
will closely match the performance (including robustness) properties of the closed-loop system with 
the centralized controller while maintaining the desired controller partitioning structure. The cost 
function is written in terms of parameters in a state- space representation of the partitioned sub- 
controllers. Analytical expressions are obtained for the gradient of this cost function with respect 
to parameters and an optimization algorithm is developed using modern computer-aided control 
design and analysis software. The capabilities of the algorithm are demonstrated by application 
to partitioned integrated flight/propulsion control design for a modern fighter aircraft in the short 
approach to landing task. The partitioning optimization is shown to lead to reduced-order subcon- 
trollers that match the closed-loop command tracking and decoupling performance achieved by a 
high-order centralized controller. 
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Introduction 


Large interconnected systems often exhibit significant coupling between the various subsystems 
thus requiring an integrated approach to controller design. Short Take-Off and Landing (STOL) 
aircraft are examples of such systems. In STOL aircraft, the forces and moments generated by 
the propulsion system provide control and maneuvering capabilities for the aircraft at low speeds 
thus creating the need for Integrated Flight/Propulsion Control (IFPC) system design. On the 
other hand, there is the need to validate certain subsystem performance and robustness character- 
istics independently of the integrated system. This calls for a level of decentralization in control 
implementation. 

One approach to integrated control design which combines aspects of centralized and decen- 
tralized control design approaches is currently being developed at NASA Lewis Research Center 
[ 1 ]. This approach consists of first designing a centralized controller, so that all subsystem in- 
terconnections are accounted for in the initial design stage, and then partitioning the centralized 
controller into separately implementable decentralized sub controllers for individual subsystems. 
By partitioning here is meant representing the high-order centralized controller with two or more 
lower order sub controllers which have input/output intercoupling such that the overall control law 
obtained on assembling the subcontrollers closely approximates the input /output behavior of the 
centralized controller. 

The objective of the current work is to describe a parameter optimization framework for obtain- 
ing partitioned subcontrollers with a decentralized hierarchical control structure. This structure is 
state-of-the-art in IFPC implementation and allows for separate implementation of airframe and 
propulsion control systems which, when assembled, will meet the overall design criteria for the 
integrated airframe/propulsion system. 

The paper is organized as follows. The controller partitioning problem with the desired par- 
titioning structure is stated first. A criterion for matching the performance of the centralized 
controller by partitioned subcontrollers while meeting the constraints on the partitioning structure 
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is then described and a “cost function” is related to this criterion. Analytical expressions for the 
derivatives of the cost function with respect to the state-space parameters of the partitioned sub- 
controllers are developed. A parameter optimization algorithm using these analytic gradients to 
minimize the cost is then described. Finally, application of the algorithm to IFPC design using 
controller partitioning is discussed for a STOL aircraft in the short approach to landing task. 


Controller Partitioning and Related Cost 


The decentralized, hierarchical controller partitioning structure is shown in Fig. 1 with ref- 
erence to an Integrated Flight/Propulsion Control system. In this figure, the subscripts and su- 
perscripts “a” and “e” refer to airframe and propulsion (engine) quantities, respectively, and the 
subscript “c” refers to commanded quantities. The interface variables z^ represent propulsion 
system quantities that affect the airframe, such as propulsive forces and moments. The structure 
is hierarchical in that the airframe (flight) controller produces commands for the interface variable 
(z ea c ) which are tracked by the propulsion subsystem. The partitioning shown in Fig. 1 is simpli- 
fied in that the controlled output errors and the interface variables z^ are assumed to be the only 
inputs to the subcontrollers (e.g. there could also be direct feedbacks for stability augmentation). 
Another simplification made in the discussion in this paper is that the plant is assumed to have no 
direct feedthrough from control inputs, i.e. the plant “D” matrix is zero. These assumptions main- 
tain simplicity while conveying the basic idea behind the formulation and solution of the controller 
partitioning problem. 

The controller partitioning problem can be stated as follows: 

Given a centralized controller with transfer matrix K(s) and a specification of the partitioning 
structure of controller inputs and outputs, i.e. 

u = K(s) e; u = 

where u a 6 R ka , u e £ R k % e a £ R m °, and e e £ R m e ; a choice of interface variables z^ £ R Pea ; a 
plant with transfer matrix G'(s) of the form 


U a 

\1 P 


, e = 




’ G(s) ' 

with 

Za 

_Gea{s) _ 


Z e 


= G(s)u and z ea = G^s) u 
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Find subcontrollers with stable transfer matrices A ra (s) and K e (s ), Le. 


u, 


= K a (s)e a and u e (.s) = K e (s) 


e c (s) 

^ea(^) 


where e ea = z ea c — Zea, so that the closed-loop performance and robustness with the sub controllers 
closely matches that with the centralized controller within the constraints of the partitioning struc- 
ture. 


The particular subsystem constraint for IFPC application is that the engine subcontroller K e (s) 
should have the structure of a command tracking controller for the interface variable commands 
z eac . Such a control structure allows the engine manufacturer to evaluate the engine subsystem 
performance independently of the airframe control and to verify that the engine subsystem will 
provide the desired performance when installed in the airframe. In general there will be practical 
constraints on the achievable bandwidth of tracking for the engine subcontroller. The lower 
bound on the Zea command tracking bandwidth will be based on achieving the desired performance 
for the integrated system, while the upper bound will be imposed by actuator limits and robustness 
requirements to high frequency modeling uncertainties. 

The approach considered in this paper for solving the controller partitioning problem is that of 
optimization of a suitable cost function over the state-space parameters of the partitioned sub con- 
trollers using an analytical expression for the gradient of the cost function. Two requirements are 
essential for the success of such an optimization approach: (1) an initial choice of parameters (par- 
titioned sub controllers) and (2) a cost function which reflects the essence of the problem statement. 
The first requirement is met via a stepwise procedure for determining partitioned sub controllers 
that meet the controller stability and the tracking structure and bandwidth requirements, as 
described in [ 2 ]. The subcontrollers obtained using this procedure are suitable for initializing the 
parameter optimization approach described herein. The cost function for controller partitioning is 
chosen to be of the form 

/(p) — «/perf(p) “I" «Arack(p) 
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where p is the vector of parameters over which optimization takes place, J pe rf(p) reflects the 
performance requirements (including robustness) and /track (p) reflects the Zga tracking requirement. 
The procedure for implementing the subcontroller stability constraint will be discussed later. In 
reference [ 3 ], a similar approach to controller partitioning was described; but, in that work, only 
the performance cost was used and a different cost function (the H°° norm) was used to express it. 

The performance cost </ pe rf(p) is chosen to be 


U?2 

Jperf = / - K(p)(ju,))G(ju;)W(ju)]) 2 au (1) 


where K(s) is the “equivalent” centralized controller obtained on assembling the partitioned sub* 
controllers using appropriate plant information (the procedure for obtaining K(s) will be discussed 
later in the paper). The above choice of / pe rf(p) corresponds to matching the loop transfer function 
matrix at the control inputs, u, during controller partitioning. Since the closed-loop performance 
and robustness are strictly a function of the loop transfer characteristics for the feedback struc- 
ture being considered here, matching the loop transfer function matrix indirectly represents the 
requirement to match the centralized controller performance and robustness with the partitioned 
sub controllers. The frequency band, [u>i,u> 2 ], is the interval over which a good match between 
loop transfer matrices is sought while the frequency dependent weighting matrix, W y allows for 
emphasizing certain frequency ranges and directions in obtaining this good approximation. In 
particular, choosing W(jcj) = (/ + /v (j^)G(ju>)) _1 , the inverse of the return difference matrix, en- 
sures that heavy emphasis will be placed on matching the closed-loop performance and robustness 
characteristics of the centralized controller by the partitioned controller. 

Note that / per f(p) is the H 2 norm of the weighted difference between the loop transfer function 
matrices. We might have preferred on theoretical grounds to use the H°° norm (cf., [ 5 ]), but this 
is numerically difficult to calculate and to differentiate in order to generate the gradient needed in 
the optimization process. 
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The Zea tracking cost /track (p) is chosen to be 

j _ 7 X' \ ( H T cent(j^) ~ TWU 

1 rack yV’l ll^centMlb 

T c *ent is the transfer function vector from the airframe commands z ac to the i th interface variable 
with the centralized controller. T* is the transfer function vector from the airframe commands 
Za c to the i th interface variable as commanded by the partitioned airframe controller, z eac l , with 
the partitioned subcontrollers and A* is a scalar weighting. || • H2 denotes the Euclidean norm 
of the row vector. Note that for the partitioned sub controllers to closely match the performance 
achieved with the centralized controller, the response of the interface variables z & to airframe 
commands z ac with the partitioned subcontrollers must match the corresponding response with 
the centralized controller because the interface variables (e.g., propulsion system generated thrust) 
significantly affect the airframe responses z a . Thus requiring T c *e n i to closely match T x , as reflected 
in /track (p)> will result in partitioned controllers such that z^ 1 appears to be a command for 
Zea* • The scalars A* are tuning parameters which allow the control designer to provide relative 
weighting for enforcing the command tracking structure among the various elements of z^ 1 and for 
trading off the performance cost against the cost of enforcing this command tracking structure. It is 
conjectured that the scalars A, provide an indirect means for adjusting the z tracking bandwidth 
achieved by the optimized engine subcontroller. For instance, choosing A* “large” will cause the 
optimization to proceed in a direction such that the z ac —► Zea* command generated by the airframe 
subcontroller “closely” matches the z ac — » z^ 1 response with the centralized controller. This 
observation together with the fact that the z a<l — ► z^ response with the partitioned sub controllers 
must closely match the corresponding response with the centralized controller in order to minimize 
/perf(p) implies that the optimized subcontrollers will be such that z^ 1 and z^* match very closely. 
This in turn implies that the optimized engine subcontroller will be such that the engine subsystem 
has “large” tracking bandwidth for the z^ 1 command. 

Cost Function and Gradient Evaluation 

Cost Function. The parameters in the optimization process are certain entries in state-space 
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realizations of K a (s ) and K e (s) as defined in the following. The notation is used throughout 
to indicate the matrix M G {A,B,C, D) in the state-space realization of the system transfer 
matrix s G {c,p,a,e} (c =centralized controller, p =plant, a =airframe subcontroller, e =engine 
subcontroller) with input i (resp. output o) G {p, a,e,ea c ,ea} 


I< a (s) : 
K'(s) : 


u e = {Ct e {sl - BU) + (Dl 



One consideration in choosing parameterizations is to introduce a “minimal” number of pa- 
rameters in the optimization process. A form for the state-space realizations introduced in [ 6 ] 
was used as the model for our parameterizations. The subcontroller system matrices A a and A e 
are represented as block diagonal matrices with two-by-two real companion blocks of the form 
. Note that, if the order of either A a or A e is odd, there is also one diagonal real entry 
corresponding to a real eigenvalue. In addition, a and 0 are constrained to be negative in order 
to meet the requirement that sub controllers be stable. In addition the first columns of each of 
the subcontroller input matrices B“ a and B e — [ B e em 5| e ] are fixed at non-zero values given by 
the initial partitioning. This last condition is different from the one used in [ 6 ] where the first 
column of the B matrix was fixed as the unit vector ej. The initial partitioning for the numerical 
example considered in this paper could not be put into this latter form necessitating the use of the 
alternative form. This alternate form is discussed in reference [ 4 ]. 

The parameters over which the optimization takes place are then the a and 0 entries in the 
block canonical forms, the entries in all but the first columns of the matrices B% a and B e and all 

c a i r D a l 

the entries in the matrices ““ , C| e , “° and [ D e eea D e ee ]. The parameter vector will be 

„'“'ea c aj ea c a . 

denoted as p G R N where N = n a (k a + m a + p^) + k a (m a + p^,) + n e (k e + m e +Pea) + (k e +Pea) m e- 
The number of parameters depends not only on the total numbers of controller inputs and outputs 
and interface variables which are fixed but also on the orders of the subcontrollers, n a and n e . 
There is thus a double incentive for keeping these orders low; not only to reduce the complexities 
of the sub controllers but also to accelerate the optimization algorithm whose performance depends 
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on the number of parameters. 


State-space representations for the transfer matrices K(s ) and (j(s) = 
while those for K(s), T cent (s) and T(s ) are constructed from the state-space representations of 
A(s), K a (s ), K e (s), Gea(s) and G(s). 


G(s) 

GM 


are given 


The transfer matrix K(s ) which enters into the performance cost term J per f(p) depends on 
K a (s ) and K e (s ), and on the transfer submatrix of the plant from control inputs (u) to interface 


u<. 

u e 


variables (z L), G m : Zea = C^ p (sl - A p ) 1 [B p a B p e ] 

assembled controller K(s ) = C(sl — A)~* B -1- D was shown in [ 3 ] to be 


A state-space realization for the 


A = 
B = 


■ (/l p - p ) ( B*fiU + B p e D 


e 

eea 


o 

ue /~ip 

eea^ eap 

’ (B p D a + B p D e D a ) 

V 1J pa^ aa ' 1J pe^eea ^ ea c a ) 

B a 

aa 

ue r\a 
D eea 1 ^ ea c a 


A a 


De rui 
eea'--' ea c a 


B p D e 

pe ee 


0 

Bl 


Qa c a) 


up fie 
u pe K - ee 


0 

A e 


C = 
D = 


— D e C p 

^ eca ' eap 

D a 

^ aa 

D e D a 

eea ea c a 


CL 

DltaC * 


Dl 



and 


The state-space representation for T cent (» = C ceni (sl- A cent ) 1 5 cent + D Qent whose rows enter 
into the evaluation of the cost term Jtrack(p) can be written in terms of the submatrices in the 


state-space representations for K (s) and G(s), i.e. K(s) 


(si - AT 1 [B c ca B c ce ) + 


and G(s ) = 


C p ap 

L C P p J 


(. si - A p ) [ B p a B p e ] . As is shown in [ 4 ], 


D c 

^ aa 

D c 

L ee 


^cent _ 
^jcent 


(-B' ca C p ap - B^Clp) 

[ (B p a CL + B p e C c ec ) A p - B p a D c aa C p v - B p e Dg e C p p 


B c 

^ ca 

B p D c 

L n pa iJ aa J 


a 


cent 


= [0 C p ) D cenl = 0. 


Similarly, T(s ) = C(sl - A) -1 B + D whose rows also appear in J t rac k(p) can be written in 
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terms of the state-space representations for K a (s), K e (s ) and G(s) 



' A a 

0 

1 

i a p 
p 

r? 

A = 

R e r a 

n eea^ea 

c a A* 

Aep 


Ap a 

B pe Ct e 

App J 



B a 

aa 



B = 


pe r\a 
u eea ^ ea c a 




(B p D a + B p D e 

Lv-^po-'^oo ' ^pe^eea 

D a ) 

^ ea c a) J 


o> 

II 

cS 

(U CO 
8 

o 

-(V%,.CS r )l and 


' ea c a 


where A ep = + B e eea C^ p + B^D^C*,), A pa = (B p a C a aa + B^D^C^) and A pp = 

(A p - B* a Dl a Cl p - B^D^D^C*, - - B p e V\ e C p v \ 

The total cost function is calculated by applying Simpson’s Rule for numerical integration over 
the frequency interval [u?i,n?2] to the sum of the two terms Jperf(p) an d Jtrack(p) described in for- 
mulas (1) and (2) respectively. The expressions developed above for the state-space representations 
of K(s),K(s),T cent (s) and T(s ) are used for calculating these costs. Note that A; in Jtrack(p) are 
to be chosen independently by the control designer to achieve the desired tracking bandwidth for 


Analytic Gradient Development. In addition to requiring the values of the “combined cost 
function,” J(p), the optimization method which is used requires the gradient V J (p) which, for any 
specific parameter vector p, is a vector with N entries consisting of the partial derivatives of J with 
respect to each parameter. The calculation of these partial derivatives for the J perf (p) portion of 
the cost function is illustrated in the following. A similar procedure is used to calculate the partial 
derivatives of /track (p)- Complete expressions for the partial derivatives are available in [ 4 ]. 

A formula for the partial derivatives of distinct non-zero singular values, 0i[A], with respect to 
an element p of the parameter vector p which can be found in reference [ 7 ] is 

= Re[u,-^ Vi l 

op dp 

where Uj and v, are right and left normalize d s ingular vectors of A and where the superscript * 
denotes the conjugate transpose. 
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For any specific parameter p (i.e. component of p), the partial derivative of Jp er f(p) with 


respect to p may be expressed as 


^^perf (p) 

dp 


d 




/ £(<Mtf( w >p)D 2dw 

CJ! 1 


dp 


U>2 

= J ^2(7vRe 


' Off( M ,p) 

U t Q V l 

dp 


du> 


Wi 


( 3 ) 


where H(u> y p) - (K(ju>) - K(p){ju))G(ju)W{ju). 

The partial derivatives of H{u, p) with respect to each of the parameters p of p are derived 
by noting that the only terms involving p in H{ u>,p) arise from the K(p)(juj) factor. Using the 
“product rule” to differentiate the state space representation, K(p)(ju) = C(ju?I — A)~ l B + -D, 
and then multiplying from the right by G(ju)W(ju>) yields 


an(u, P ) _ ac( p) WuJ _ X) . , SGUw)w{ju) 


dp 


dp 


- C(ju>I - ^l(p))- 1 - A{p))- l BGVu)W(ju) 

- C (ju>I - A)- 1 ^M G {ju)W{ju) - ^ '±G(ju)W(ju > ). (4) 


Substituting (4) in (3) yields an expression for 


dp 

9 </perf (p) 

dp ‘ 


dp 


The final step needed to carry out the determination of the gradient of J pe rf is to calculate 
the partial derivatives of A, 5, C, and D with respect to the parameters which appear in these 
matrices via the sub controller state- space matrices 


M <E 


as described earlier. Each parameter p corresponds to a specific row and column element of some 
parameter submatrix, say the element in row j y column k of submatrix M. Any block of A, B, etc. 
not containing M as a factor will be replaced by a zero block in the calculation of the corresponding 
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partials — , etc. On the other hand, any block containing M as a factor contributes the same 

dp dp 

block with the factor M replaced by the submatrix e_,e fc T , which has the j,k th - entry as its only 
nonzero element. Here Gj is the vector whose j ^ element is 1 and all other elements are zero. For 
example if p corresponds to the row j, column k entry of A a , then 

0 0 O' 

0 ejGk T 0 
0 0 0 


8A 

dp 


hile - 7 T- and are all 0. If p corresponds to the row j, column k entry of Dl ea , then 
dp dp dp 


dAjp) 

dp 


0 0 0 
0 0 0 


with expressions for and derived in a similar manner. Here the 0 entries 

dp dp dp 

refer to zero blocks. These two examples demonstrate extremes of complexity in the calculation of 
the various partials of the state space matrices for K{s). Once these partials have been determined 
with respect to all parameters in p, the resulting matrices can be inserted into the formula for the 
partials of J per f described earlier and the gradient can be calculated by applying Simpson’s Rule of 
numerical integration over the same frequency interval as was used for evaluating J peT t(p)- 

The partial derivatives of Jtrack are computed by a similar procedure and these are combined 
with the partials of J per f(p) to yield VJ(p) = VJ per f(p) + V«7 t rack(p)- 


The Optimization Algorithm 


The flow of the optimization algorithm for controller partitioning is shown in Fig. 2. 

The fixed data used by the algorithm are state-space representations for the plant submatrices 
G(s) and G ! ea( 5 ) f the centralized controller K(s), the weighting matrix VF(s), as well as a parti- 
tioning structure for the controller inputs, outputs and interface variables. The control designer 
must also input values of the tracking weight parameters A; used to determine T tC ack(p)- The main 
steps in the algorithm are: 
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1. The initial partitioning ( K$(s ) and Kq(s)) is obtained by applying the stepwise procedure 
described in [ 2 ]. Special attention is paid to obtaining reasonably low-order sub controllers 
which are stable and satisfy the command-tracking requirement. 

2. The initial partitioning is converted to the “minimal parameter” form to generate an initial 
value of the parameter vector. 

3. The (initial as well as any subsequent value of the) parameter vector is passed to a module 
which determines the state-space representation for the hierarchically partitioned controller 
and calculates the combined cost. The gradient is also computed analytically by the procedure 
described above. 

4. The current gradient is used in conjunction with previous information to generate a direction 
of search. A one-dimensional linesearch is carried out using the cost and gradient calculated 
at each parameter point to predict a new parameter point until one is found which yields a 
sufficient reduction in both the cost function and the size of the gradient. This linesearch is 
constrained so as to maintain stability of the subcontrollers. 

5. At the end of the linesearch, the new parameter point, cost and gradient are compared to the 
values at the beginning of the linesearch as a check on convergence. If either the maximum 
absolute value of the partial derivatives is less than a user-specified tolerance or both the 
maximum change in all the parameters is less than one user-specified value and the change in 
the total cost is less than another value then convergence is declared and the iteration ceases. If 
the convergence test fails, the algorithm proceeds to update the information used to determine 
the direction of search and to use the most recent cost and gradient values to generate a new 
direction of search and carry out the linesearch via steps 3. and 4. 

6. The output of the algorithm is the state-space representation for sub controllers which yield the 
assembled hierarchically partitioned controller with performance and stability robustness com- 
parable (within convergence criteria) to that of the centralized controller and with reasonable 
adherence to the command-tracking requirement. 
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7. These subcontroller transfer matrices have the same orders n a and n e as the initial partitioning. 

Controller reduction can be performed on these “optimal subcontrollers” and the process of 

optimization can be repeated on the “new initial partitioning”. 

We use the Broyden Fletcher Goldfarb Shanno implementation of the quasi-Newton method 
to carry out this optimization, see [ 8 ]. This iterative method requires the calculation of the 
combined cost and gradient for parameters as referred to above. It uses successive gradients to 
build up an approximation to the inverse Hessian matrix. In this way, subsequent search steps are 
successively closer to those generated by Newton’s Method and convergence is accelerated as the 
iterations proceed. 

The results of applying this process to an IFPC example are reported in the following section. 

Example of Controller Partitioning 

The controller partitioning algorithm was applied to the centralized flight /propulsion controller 
for a STOL aircraft obtained in reference [ 9 ]. This controller has the form u = K(s)e with the error 
vector e consisting of errors, e = [e q , e„, fijV2, Cepr] T , in following pitch rate variable (q v = <7+0.1$), 
velocity (v), engine fan speed (N 2) and engine pressure ratio ( EPR ) commands. The control input 
vector u consists of rates of change of thrust vectoring angle, fuel flow, thrust reverser port area 
and nozzle throat area, u = [6 tv> WF, A78, A8]. Note that u consists of rates because integrators 
were appended to the control inputs during the process of centralized control design to achieve 
zero steady-state error for step commands. Based on open loop control effectiveness studies for the 
plant, the partitioned airframe and engine controllers are desired to have inputs e a = [e q ,ev] T and 
e e = [cat 2 > C£ph] t and outputs u tt = [$xv] a nd u e = [ WF, A78,j48] t respectively. The integrated 
plant and the centralized controller state-space matrices are listed in the Appendix. The interface 
variable for this design example is FEX, the axial thrust generated by the propulsion system, 
since the main function of the propulsion system is to provide the necessary axial thrust required 
to track pitch rate and velocity commands. The structure of this controller partitioning is shown 
in Fig. 3. 
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An initial controller partitioning was obtained using the procedure discussed in [ 2 ]. The plant, 
centralized controller, and initial partitioning state-space matrices are listed in the Appendix. The 
airframe and engine sub controllers for this initial partitioning are of order 10 and 7, respectively. 
The performance of the initial controller partitioning was evaluated in comparison with that of the 
centralized controller by comparing closed-loop system response to step commands in the controlled 
variables z. The responses to q Vc , N 2 C and EPR C with the initial partitioned sub controllers were 
comparable to those with the centralized controller. However, the response to V cy shown in Fig. 
4, showed considerable degradation in terms of increased coupling in the N 2 and EPR responses 
with the partitioned subcontrollers. Note that all the quantities shown in Fig. 4 were normalized, 
using scalings discussed in [ 9 ], to allow a direct comparison of the various response magnitudes. 

For all the results presented in this paper, the optimization was done over the frequency range 
u) E [0.1,100] with 20 frequency points per decade. The number of parameters in the optimization 
was 95 with 44 corresponding to the airframe subcontroller and 51 corresponding to the engine 
sub controller. The algorithm was first applied with A = 0 ( p ea = 1, so A; = A), i.e. J(p) — J per f(p) 
with no attention to the z e0c command tracking structure requirement for the engine subcontroller. 
As seen from Fig. 5, the optimization code showed good numerical behavior with convergence to 
a significantly reduced cost J(p) in 100 iterations. The maximum singular value of the weighted 
loop transfer error — K )G{I + KG)“ l ](ju>), which is a major contributor to the cost 

J(p) in this case, shows a considerable improvement with partitioning optimization as is seen from 
Fig. 6. The results in Fig. 6 imply that the optimized partitioned controllers will more closely 
match the performance and robustness characteristics of the centralized controller than will the 
initial partitioned controllers. The closed-loop response to step V c with the optimized partitioned 
controllers is shown in Fig. 4 and is comparable to that for the centralized controller both in terms 
of tracking of the velocity command and decoupling of q v , N 2 and EPR responses. The closed 
loop responses with the optimized partitioned controllers for step g„ c , N 2 C and EPR C were also 
checked to ensure that these were comparable to the responses with the centralized controller. 

Figure 7 shows the Bode gain plots for the interface variable command tracking response 
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( \FEX/FEX C \ ) for the engine subsystem with the initial partitioning and the optimized partition- 
ing obtained using A = 0. The engine subcontroller for this optimized partitioning has a FEX 
command tracking bandwidth of approximately 4 rad/s (tracking bandwidth is defined as the fre- 
quency at which the Bode gain is —3 dB), which is quite low and implies that the engine subsystem 
will not provide adequate tracking of FEX C in the presence of disturbances and low frequency 
variations in the engine dynamics. In order to verify command tracking and disturbance rejection 
properties of the engine subsystem in independent subsystem check-out, it is desirable to have a 
minimum bandwidth of around 10 rad/s for FEX command tracking. Since FEX has significant 
control effect on the aircraft responses V and q v , such a lower bound on FEX command tracking 
is also required to guarantee performance robustness of the integrated system in terms of tracking 
and decoupling of responses to V and q u commands. As seen from Fig. 7, the initial controller 
partitioning satisfies this FEX command tracking bandwidth. 

The fact that the optimized partitioning with A = 0 gives performance that is comparable to 
the performance achieved with the centralized controller although the optimized engine subcon- 
troller has poor FEX command tracking properties can be understood by analysis of the results 
shown in Fig. 8. We note that the FEX response to V c with the optimized (A = 0) partitioned 
sub controllers closely matches the corresponding response with the centralized controller. This 
result is as expected because FEX significantly affects the V and q v response of the aircraft and 
matching the FEX response will be required if the partitioned subcontrollers are to match the 
performance achieved with the centralized controller. However, as seen from Fig. 8, the command 
for the interface variable FEX C as generated by the partitioned airframe subcontroller significantly 
“leads” the actual FEX response. So, the optimized airframe subcontroller compensates for the 
low interface variable command tracking bandwidth property of the engine subcontroller by gener- 
ating a command for the interface variable which leads the desired interface variable response by an 
appropriate amount. This phenomenon of lead compensation for low tracking bandwidth is further 
evident from the Bode gain plots shown in Fig. 9 for the FEX/V C response with the centralized 
controller and the FEX/FEX C , FEXJV and FEX/V C responses with the optimized controller 
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partitioning (A — 0). In Fig. 9, the Bode gains of FEX C /V C ^nd FEX/FEX C do not exactly 
add up to the Bode gains of FEX/V C for the partitioned subcontrollers because the FEX/FEX C 
response is only for the engine subsystem and does not account for the effect of thrust vectoring 
(u a ) on the FEX response. 

As previously discussed, there is a minimum z^ command tracking bandwidth which must be 
provided by the partitioned engine subcontroller for independent validation of the engine subsystem 
and performance robustness of the integrated system, and there is an upper limit on the achievable 
bandwidth imposed by actuator and stability robustness constraints. Some “lead in the interface 
variable command Zea c generated by the airframe subcontroller is desirable to compensate for the 
limited bandwidth engine subsystem. Nevertheless, large “lead” is undesirable as it can result in 
quick saturation of the engine actuators because of command magnification. The effect of varying 
the user specified weights (A; in equation (2)) in the controller partitioning optimization procedure 
to address the trade-off between lead and bandwidth is demonstrated by the plots of the Bode 
gains of the FEX / FEX c responses shown in Fig. 7. As is seen there, the FEX command tracking 
bandwidth achieved by the optimized engine subcontroller increased as A was increased. The closed- 
loop system performance of the optimized partitioning with A = 0.05 for step V c closely matched 
the performance with the centralized controller as shown in Fig. 10 (note that the scales in Fig. 10 
are different from those in Fig. 4). The airframe subcontroller obtained from the A = 0.05 optimal 
partitioning also provided small lead in the FEX command, as is seen from Fig. 11, and the FEX 
command tracking bandwidth of 20 rad/s (from Fig. 7) for the corresponding engine subcontroller 
was within the engine actuation limits. The state-space matrices for these airframe and engine 
sub controllers are listed in the Appendix. 

With the optimized controller partitioning for A = 0.05 as a baseline, the possibility of reducing 
the order of the partitioned subcontrollers was investigated. The engine subcontroller was reduced 
to order by residualization of the three high frequency modes without any loss of performance. 
Through the use of internally balanced reduction techniques [ 10 ], the airframe subcontroller was 
reduced to 6 th order (from the original 10 th order) without excessive mismatch in the controller 
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transfer matrix characteristics as is seen from the full and reduced order airframe controller singular 
values comparison in Fig. 12. This reduced order airframe subcontroller did, however, exhibit de- 
terioration in closed-loop performance as shown in the V and q v response comparison plots for step 
V c in Fig. 13. The controller partitioning optimization procedure was performed with A = 0.05 and 
with optimization over only the airframe subcontroller parameters (the engine subcontroller was 
fixed to the reduced 4 th order controller). The response obtained with the optimized reduced order 
airframe subcontroller for step V c is also shown in Fig. 13. Note that the optimized subcontroller 
provides improved tracking of the velocity command. Although, as can be seen from Fig. 13, the 
transient in q v is larger in magnitude for the optimized reduced order subcontroller than for the 
initial reduced order airframe sub controller, the transient is of significantly shorter duration and 
the pitch variable quickly settles down to the zero commanded value. The state-space matrices 
for the optimized reduced order subcontrollers are listed in the Appendix. These results demon- 
strate the potential for using the controller partitioning procedure in conjunction with controller 
order reduction techniques to obtain the simplest (lowest order) partitioned sub controllers that will 
provide adequate integrated system performance. 

The results presented so far have focused on comparing the performance achieved with the op- 
timized partitioned sub controllers with that achieved with the centralized controller. Robustness 
issues are also of importance in practical control design. Robustness analysis was performed using 
structured singular values for gain and phase variations occuring in the controlled outputs. The 
results are shown in Fig. 14 for the centralized controller, the optimized full order partitioned sub- 
controllers corresponding to A - 0.05 and the optimized reduced order subcontrollers.The procedure 
for creating the interconnection matrix to perform gain and phase margin robustness analysis using 
structured singular values is documented in Ref. [ 11 ] and other references listed therein. From 
Fig. 14 the stability margin parameter /z corresponding to the maximum over frequency of the 
structured singular values is 1.37 for the centralized controller and 1.43 for the full and reduced or- 
der optimized partitioned subcontrollers. These values of /i translate into guaranteed multivariable 
gain margins of —4.8 dB to 11.4 dB and —4.6 dB to 10.5 for the closed-loop system with the cen- 
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tralized and partitioned sub controllers respectively, and similarly guaranteed multivariable phase 
margins of ±42.8 and ±41.0 deg, respectively, for simultaneous gain or phase variations occuring 
in all loops at the controlled outputs. These results indicate that the robustness characteristics of 
the centralized controller are maintained by application of the parameter optimization approach. 

Conclusions 

A parameter optimization framework was presented to solve the problem of partitioning a 
centralized controller into a decentralized, hierarchical structure suitable for integrated flight/ 
propulsion control implementation. A cost function to be minimized was formulated such that 
the “equivalent” centralized controller assembled from the “optimized” partitioned subsystem con- 
trollers will closely match the performance and robustness of the closed-loop system with the 
centralized controller while maintaining the desired characteristics of the controller partitioning 
structure. An analytical gradient based optimization algorithm was developed and applied to con- 
troller partitioning of an integrated flight /propulsion control design for a modern fighter aircraft 
in the longitudinal short approach to landing task. The optimization of the partitioned subcon- 
troller state-space parameters resulted in improved performance over an initial partitioning and 
close matching of the performance with the centralized controller. The ability of the optimization 
to address the subsystem control structure requirement of an adequate command tracking band- 
width for the propulsion subsystem by varying scalar weight parameters in the cost function was 
demonstrated for the example study. It was shown that by judicious use of the partitioning opti- 
mization procedure in conjunction with model reduction techniques, simplified low-order subsystem 
controllers with a desired structure can be obtained that match the performance and robustness 
characteristics of a high-order centralized controller. 
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Appendix 


The data matrices and initial partitioning matrices for the controller partitioning example are listed below. In all cases 
the given matrices A, B, C and D correspond to the state-space representation of the given system or subsystem 

dx , „ 

ii = Ax+B>r 


u = Cx + Dy. 


These matrices will be given in the form 


A 

B' 

C 

D 


if space allows or will be written out individually if the former is too large. 

The integrated airframe propulsion system with integrator augmentation is represented by 
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The centralized controller has state-space representation 
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The initial partitioning consisted of and airframe controller with representation SKao of order 10 and an engine controller 
with representation SKeo of order 7. The initial partitioning was obtained by the procedure of reference [ 2 ]. This partitioning 
is not in the minimal parameter form. These are split because of their size. 
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The initial engine controller is represented by 
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The parameter optimization algorithm for controller partitioning was applied to the problem with the starting matrices 
SKao and SKgo listed above. The controllers obtained from this process had state-space representations SKAopt for the 
airframe and SKEopt for the engine. 
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The optimized engine controller has state-space representation 
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Residualization of high frequency modes was applied to the optimized engine controller to reduce it to one with order 4. 
Balanced model reduction was applied to the optimized airframe controller to reduce it to one of order 6. The optimization 
procedure was applied to this sixth order subcontroller with the engine controller fixed at the one of fourth order. The 
resulting reduced order optimized airframe subcontroller is 
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Figure 1 . — Controller partitioning structure. 


Figure 2. — Flowchart for partitioning optimization alagorithm. 
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Figure 3.— Partitioned controller structure for design example. 





Figure 4. — Closed-loop system response to step velocity 
command for centralized controller, initial partitioning and 
optimized partitioning with X = 0. 
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igure 5. — Cost 
with X = 0. 


for controller partitioning optimization 
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Figure 6.— Maximum singular value of weighted error, 
a maxt(K-K)G(l+KG) _1 ]{jw) p for initial and optimized (X = 0) 
partitionings. 
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Figure 7. — Bode gain plot of engine closed-loop system 
response (KEXc — FEX) for initial and optimized partitionings 
with varying X. 



Figure 8. — FEX response to step velocity command for central- 
ized and optimized partitioning with X = 0 and FEX** response 
for optimized controller (X = 0). 
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Figure 9. — Lead vs. bandwidth trade-off in controller 
partitioning. 




Figure 1 0.— Closed-loop system response to step velocity 
command for centralized controller and optimized partitioning 
with X = 0.05. 
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Figure 1 1 . — FEX response to step velocity command for 
centralized controller and optimized partitioning with X = 0.05 
and FEXc response for optimized controller (X » 0.05). 
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Figure 13. — Closed-loop system response to step velocity 
command for centralized controller and reduced order 
partition subcontrollers { n a = 6 n e = 4) - initial and opti- 
mized (K e fixed). 
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Figure 1 2. — Singular values of the airframe subcontroller for 
controller partitioning with X = 0.05 - full (1 0) and reduced (6) 
orders. 
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Figure 1 4. — Structured singular value for robustness analysis 
to gain and phase variations at the controlled outputs z with 
centralized controller, for optimized partitioned controller 
(X = 0.05), and optimized (K® fixed) reduced order controller. 
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